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ABSTRACT 


An  automatic  station -keeping  satellite  must  be  able  to  measure 
its  deviation  from  the  nominal  position.  For  the  system  considered 
here  (an  equatorial  synchronous  orbit)  this  involves  a  latitude  differ¬ 
ence  measurement  between  the  earth  and  sun  at  certain  times.  This 
report  shows  that  a  combination  of  an  error  in  satellite  orientation 
and  the  sun  being  in  a  plane  other  than  the  equatorial  plane  may  re¬ 
sult  in  errors  in  position  determination.  The  nature  of  the  errors 
involved  is  described  and  their  magnitudes  estimated. 
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Position  Error  in  Station-Keeping  Satellite 


Ideal  System  Operation 

The  satellite  is  in  a  synchronous  equatorial  orbit  and  its  spin  axis  is  par¬ 
allel  to  the  orbit  normal.  It  is  equipped  with  sensors  which  measure  the  lon¬ 
gitude  difference  in  satellite  coordinates  between  the  sun  and  the  earth.  If  the 
satellite  is  accurately  synchronous,  the  longitude  difference  should  be  a  cer¬ 
tain  value  (say  90°)  at  the  same  time  (measured  by  an  on-board  clock)  each  day. 
Thus  a  measurement  of  the  time  at  which  the  longitude  difference  is  90°  may  be 
used  to  control  a  position  correcting  mechanism. 


Statement  of  Problem 

How  will  this  system  function  when  the  satellite  spin  axis  is  not  parallel 
to  the  orbit  normal  and  the  sun  is  not  in  the  plane  of  the  equator?  That  is, 
what  errors  may  be  expected  in  the  estimation  of  longitude  difference? 

I  first  attempt  rough  estimates  of  the  errors  involved  on  an  intuitive  basis. 
A  computer  program,  based  on  the  vector  approach,  is  described  and  results 
from  this  substantiate  the  numbers  of  the  next  section. 
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Summary  of  this  Note 


Due  to  the  sun's  elevation  an  error  9  in  satellite  spin  axis  position  may 
(depending  on  the  details  of  the  geometry)  couple  into  an  error  in  orbit  posi¬ 
tion  estimation  of  9  x  .  4.  It  appears  at  this  date  that  an  error  of  this  magni¬ 
tude  is  quite  significant  to  functioning  of  the  position  control  mechanism. 

This  may  indicate  that  a  reduction  of  attitude  error  is  required. 

Geometry  of  the  Problem 

1.  First  consider  the  problem  when  the  sun  S  is  in  the  equatorial  plane  of 
the  earth  as  in  Fig.  1  and  the  satellite  axis  A*  makes  an  angle  0  with  the  orbit 
normal.  Now  E*  S'  is  the  longitude  difference  measured  by  the  satellite.  If 
the  satellite  were  correctly  oriented  at  this  point  in  the  orbit,  it  would  meas¬ 
ure  the  arc  ES.  Let  OL  be  the  line  of  nodes.  Then  from  spherical  geometry 


tan  LE 


tan  LE1 

cos  9 


tan  LS 


tan  LSf 

cos  9 


ES 


tan 


-1 


f  tan  LS,n\ 
x  cos  0  J 


+ 


tan 


•  1  tan  LE'  A 

^  cos  9  J 


Thus  for  0  <  10°  the  error  is  less  than  l1. 

2.  Now  consider  the  effect  of  the  sun's  declination  on  the  problem.  The  pre¬ 
vious  paragraph  has  shown  that  with  the  sun  in  the  earth's  equatorial  plane  the 
satellite  measures  the  arc  E'S'  —  ES.  Now  let  the  sun  be  at  its  zenith  Sj  .  Let 
the  spin  axis  move  from  N  to  A*:  then  for  small  tilt  angles  we  may  regard  the 
sun  longitude  line  A'S^S^'  as  hinging  at  S^.  In  order  to  estimate  this  lever  arm 
effect  we  project  onto  a  plane  perpendicular  to  the  line  OS^.  By  the  arguments 
of  the  previous  paragraph  the  error  in  longitude  estimation  is  SS^~S^S'.  We 
find  SS^  —.4x0  for  the  worst  case.  The  worst  case  is  when  the  right  ascen¬ 
sions  of  the  sun  and  spin  axis  differ  by  odd  multiples  of  90°.  The  error  ap¬ 
proaches  zero  when  the  right  ascensions  differ  by  even  multiples  of  90°, 
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Computer  Program  Results 


These  are  plotted  on  Figs.  Z,  3,  and  4.  On  each  graph  there  are  Z  curves. 
One  curve  is  for  the  sun  at  zenith  and  the  other  for  the  sun  is  the  earth's  equa¬ 
tor.  The  position  of  the  sun  is  fixed  for  each  curve.  The  longitude  difference 
between  earth  and  sun  as  seen  by  a  correctly  oriented  satellite  at  that  point  at 
which  a  tilted  satellite  measures  this  difference  as  90°  is  plotted  as  ordinate. 
The  curve  is  now  obtained  by  letting  the  spin  axis  move  around  the  North 
Pole  at  a  constant  angle  9  and  calculating  this  longitude  difference  at  each  point. 
The  graphs  are  for  0  =  Z0°,  1°,  3°  respectively. 

Vector  Solution  of  the  Problem 


3-63-4747 


VERNAL 

EQUINOX 


EARTH 


Given:  The  right  ascensions  and  declinations  of  the  spin  axis  A,  the  sun  S, 
the  longitude  difference  in  satellite  coordinates  between  the  sun  and  earth. 

Required:  To  find  the  longitude  difference  that  would  be  measured  by  a  cor¬ 
rectly  oriented  satellite  at  this  point. 

Compute  coordinates  of  unit  vectors  A,  S  from  the  right  ascensions  and  dec¬ 
linations  . 

The  X,  Y,  Z  components  of  a  unit  vector  V  may  be  obtained  from  the 


right 

ascension  and 

vi 

=  cos  fi 

cos  X 

V2 

=  cos  0 

sin  X 

V, 

=  sin  0 

3 
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We  obtain  components  for  the  two  vectors  S  and  A. 

We  now  obtain  the  position  of  E  given  the  positions  of  S  and  A  and  that  the 
longitude  differences  in  satellite  coordinates  between  S  and  E  is  cj). 

The  vector  A  x  E  is  in  the  equatorial  plane  of  the  satellite,  at  right  angles 

to  the  plane  containing  A  and  E.  Let  AE  be  the  angle  between  A  and  E.  Thus 

A  x  E  =  (j,  _  sin  AE  where  n  is  a  unit  vector. 

E  E 

Similarly  AxS  =  fx  ^  |  A  x  S  |  . 

Now  a  positive  rotation  <j>  about  A  will  bring  fi  into  coincidence  with  fj,  . 

/i  iL  u 

M  E  X  n  s  =  sin  (j)A  A 


We  also  have 


But 


Ai  E  •  M  s  =  cos  ^ 


HS  x  (ME  xms)  =  HE  -  (ME  •  MS)MS 


or 


M  E  =  M  s  x  A  sin  (j>A  +  cos  (}>A  n 


A  S 


This  equation  may  be  solved  to  give  the  components  of  (X^.  E  is  the  vec- 

E 

tor  which  fulfills  the  following  conditions. 

1.  E  •  me  =  0 

2.  E^  =  0  (the  satellite  is  in  an  equatorial  orbit) 

3.  | E  |  =  1 

4.  A  x  E  =  iu  E  sin  AE  when  AE  <  180° 
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We  obtain 


E 


1 


M 


E2 


+  /(MEi)Z  +  (Mr,)2 


El 


E 


2 


(M  ei> 


We  now  calculate  the  right  ascension  X  of  the  earth. 

E 


cos  X £  =  E^ 


sin  =  E 


tan 


2 

sin  X £ 

1  +  cos  XT 


2 

1  T  E. 


Let  be  the  right  ascension  of  the  sun.  Then  at  this  point  in  the  orbit,  when 

the  tilted  satellite  measured  a  longitude  difference  between  the  earth  and  the 

sun  of  90°  a  correctly  oriented  satellite  would  measure  -  X 

b  E 
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Fig.  1.  Correctly  oriented  satellite  (spin  axis  at  A)  measures 

arc  ES  (-  90°).  With  sun  at  S  the  satellite  measures  E'S'  ~  ES. 

With  sun  at  S.  the  satellite  measures  E'S!,  ~  ES  -  .4x9. 

1  1 
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NORTH  POLE 


(&ap)  aoaa3  Nomsod 
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RIGHT  ASCENSION  OF  SPIN  AXIS  (deg) 


NORTH  POLE 


(Bap)  dOaa3  NOIllSOd 
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RIGHT  ASCENSION  OF  SPIN  AXIS  (deq) 

Fig.  3 


NORTH  POLE 


(  6»P  )  HOHH3  NOIllSOd 
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RIGHT  ASCENSION  OF  SPIN  AXIS  (deq) 

Fig.  4 


APPENDIX  I 


Computation  of  Position  Error 

1 .  Input 

(a)  Right  ascension  and  declination  of  the  spin  axis  A(RA,  DA) 

(b)  Right  ascension  and  declination  of  sun  S(RS,  DS) 

(c)  Earth-sun  longitude  difference  L.ES  PRINT  RA,  DA,  RS,  DS 

2.  Compute  components  of  A  along  the  axis 

A1  =  cos  DA*  cos  RA 

A 2  =  cos  DA*  sin  RA 

A3  =  sin  DA 


3.  Compute  components  of  S  along  the  axes 

51  =  cos  DS*  cos  RS 

52  =  cos  DS*  sin  RS 

53  =  sin  DS 

PRINT  Al,  A2,  A3,  SI,  S2,  S3 

4.  Compute  vector  A  x  S 
AS1  =  A2*S3  -  S2*A3 
AS2  =  A3*S1  -  S3*A  1 
AS3  =  Al  *S2  -  S 1  *  A  2 


5.  Compute  unit  vector  US 

Form  SQAS  =  (  (AS1)2  +  (AS2)2  +  (AS3)2  J1/2 


US1 


AS1 

SQAS 

AS2 


USZ 


SQAS 


US3 


AS3 

SQAS 


V 


PRINT  AS1 ,  AS2,  AS3 
US1,  US2,  US3 

Calculate  MAGUS  =  (  (US1)2  +  (US2)2  +  (US3)2)*^2  and  PRINT  MAGUS 

6.  Compute  unit  vector  UE 

UE1  =  (US2*A3  -  A2*US3)*  sin  LES  +  US1*  cos  LES 

UE2  =  (US3*A1  -  A3*US1)*  sin  LES  +  US2*  cos  LES 

UE3  =  (U S 1  *A  1  -  A1*US2)*  sin  LES  +  US3*  cos  LES 

PRINT  UE1,  UE2,  UE3 

and  ^(UEl)2  +  (UE2)2  +  UE3)2 


7.  Compute  vector  E 

+  UE2 


El  = 


+  /  (UE2)2  +  (UE1)2 


A. 


A. 


E2  = 


-  UE  1 


+  '  (UE2)2  +  (UE  1  )2 


A. 


A. 


- 2 - T 

PRINT  El,  E2  and  y  (UE2)  +  (UE1) 


8.  Compute  right  ascension  of  the  earth 

RE  =  2  tan  *  +  27T  ,  mod  27 r 

convert  to  degrees  and  PRINT 

9.  Compute 

X  =  RS  -  RE  +  360,  mod  360 
PRINT  R.  A.  Diff.  =  X 
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